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Sensitivity  of  MINQUE  with 
Respect  to  A Priori  Weights 


of  a priori  weights  needed  for  the  estimators.  The  MINQUE  have  minimum 
variance  in  the  class  of  unbiased,  quadratic  estimators  when  (i) 
and  are  normally  distributed  and  (ii)  the  ratio  of  the  a priori 
weights  is  equal  to  the  ratio  of  the  variances  of  or  and  e^.  This 
paper  shows  that  for  the  a priori  weights  in  a neighborhood  of  the 
ratio  of  the  variance  components  the  variances  of  the  MINQUE  are  quite 
insensitive.  That  is,  the  variances  deviate  little  from  the  optimum 
variances. 

For  comparison,  the  variances  of  the  corresponding  Henderson 
Method  I or  analysis  of  variance  type  estimators  of  the  variance  com- 
ponents are  also  computed.  Recommendations  can  be  made  as  to  when  to 
use  which  type  of  estimator  and  how  to  specify  the  a priori  weights  if 
a MINQU  estimator  is  used. 
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2 2 

of  o and  o is  given  as  the  solution  of  So  * u (Rao  [1972])  where  the  matrix 
S and  the  vector  u will  be  defined  below.  Both  S and  u depend  on  the  matrix 


, V*_1J  v*-1 

V*"1  - 2 

(l'V*_1l) 


where  V* 


w V,  + w 1,  J is  an  (n  x n)  matrix  of  ones  and.  w and  w are 
a 1 e n a e 


a priori  weights  which  are  needed  for  the  estimator  (actually,  the  estimator 

w 

depends  on  the  a priori  weights  only  through  the  ratio  — ) . 


1 


2 


The  problem  considered  here  is  that  when  a and  c are  normally  dis- 
tributed the  MINQU  estimators  have  minimum  variance  in  the  class  of  un- 
biased quadratic  estimators  only  when  the  ratio  of  the  a priori  weights 

w o2 

equals  the  ratio  of  the  variance  components,  i.e.  — = . At  first 

w 2 

e °e 

glance,  it  would  seem  that  this  is  a very  severe  requirement  to  assure 

a desired  optimal  property  of  the  unbia-%d  quadratic  estimator;  however, 

'2  A 2 

it  was  found  that  the  variances  of  a and  o remained  close  to  the  mini- 

a e 9 

w az 

mum  variances  for  — in  a large  neighborhood  of  — — for  several  values  of 
W 2 

e °e 

the  parameters  and  several  degrees  of  unbalancedness  of  the  design.  This 

is  discussed  in  detail  in  section  III. 

For  comparison,  the  variances  of  the  Henderson  Method  I estimators 
2 2 

of  o and  o were  computed.  These  quadratic,  unbiased  estimators  are 

6 d 

easily  calculated  and  it  is  of  interest  to  note  the  magnitude  of  their 

variances  relative  to  the  minimum  variances  attainable  in  the  class  of 

unbiased,  quadratic  estimators  as  well  as  the  variances  of  the  MINQU 

w 

estimators  which  depend  on  the  ratio  — . 

w 

e 2 2 

The  Henderson  Method  I estimators  of  a and  o are  given  as  (Searle 

a e 


[1968]) 


'2  _ n[ (n-i)«SSA  - U-1)«SSE] 

a 2*2 
(n-iHn  - l r ) 

i-1  1 


"2  SSE 
a m — ■■  - ■■ 
e (n-1) 


respectively  where  SSA  and  SSE  will  be  defined  below. 


It  is  iiqportant  to  note  that  SSA  and  SSE,  and  hence  the  variances  of 
*2  *2 

o and  o , do  not  depend  on  a priori  weights  as  do  the  MINQU  estimators. 

A A 
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IX.  Variances 

For  the  random  effects  model  described  above,  the  S matrix,  u 
SSA,  and  SSE  are  defined  as  follows: 

s - lstA 


2*2 


ij1 


where  Sn  - trtRV^V^,  S12  = S21  = trOtt^R),  and  S22  = tr  (RR) 

. - m 

2*1  LU2 J 

where  = Y'RVjRY  and  u2  = Y'RRY, 

SSA  * Y'  fU.  (U.'U,  )“ V'--  J ]Y  = Y'Q  Y, 

- Ill  Inn-  - a- 


and 


SSE  = Y'  [I-U,  (U'U,  ) U,'  ] Y = Y'O  Y. 

— ill  1 * “ e- 

2 2 

The  MIN9U  estimators  of  a and  a are  then  given  by 

a e 


and 


i2  S12U2~S22U1 

3 " S2  -S  S 
S12  11  22 


'2  „ S12U1~S11U2 

e s2  -s  s 

12  11  22 


respectively. 


*2  *2 

The  variances  of  o and  a can  now  be  found  as  follows: 
a e 

*2  2 2 2 

var(Oa)  - k {S^var  (u2>+S22var  (u^^S^S^cov  (u^Uj)  } 


and 


*2  2 2 2 

var <©e)  - k {s12var(u1)+S11var(u2)-2S12S11cov(u1,u2>} 

2 -1 

where  k ■ *S12~S11S22^  * U8*n9  results  as  found  in  Searle  ([1971] 


vector , 


pg.  55) 


and  noting  that  l'R  = O',  the  following  quantities  can  be  obtained: 
var (u^)  - 2tr (RV1RVRV1RV) , 
var(u2>  = 2tr (RRVRRV)  and 
covfu^u^  = 2tr(RVxRVRRV) 

where  V is  the  variance  of  Y as  defined  above.  Hence,  the  variances  of 

the  MINQU  estimators  can  be  written 

*2  2 2 2 2 2 

var  (°a>  = 2k  {S12tr[(RRV)  l+S^tr  [ (RV^V)  l_2S22Si2tr fRVlRVRRVJ * 

and 

var(o^)  53  2k2{S^2tr[(RV1RV)2]-fS^1tr[(RRV)2)-2S12S11tr[RV1RVRRV)}. 
"2  "2 

The  variances  of  oa  and  o^,  the  Henderson  Method  I estimators,  can 

be  found  as  above.  Noting  that  the  quadratic  forms  SSA  and  SSE  are  in- 

~2  "2 

dependent,  the  variances  of  c and  a can  be  written  as  follows: 

a e 

- 2 n2  2 2 

var (o  ) - { (n-i)  var (SSA)  + (1-1)  var(SSE)} 

/ „>2,  2 r 2»2 

(n-i)  (n  - E r. ) 

i=l  1 

and 


"2  1 

var(o  ) ■ ■?  {var (SSE)}. 

e , .,2 

(n-l) 

Again  using  results  as  found  in  Searle  [(1971)  pg.  55)  and  noting  that 

l'Q  - l'Q  ■ o',  then 
- a - e - 

var (SSA)  - 2tr(Q  VQ  V) 

A A 

and 


var  (SSE)  - 2tr(Q  VQ  V) . 

e e 


Therefore,  the  variances  of  the  Henderson  Method  I estimators  can  be  written 
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var (a  ) = 
a 


2n 


, . , 2 . 2 : 2 2 
(n-£)  (n  - £ r . ) 

i-1  1 


{ (n-£) 2tr (Q  VQ  V)+(£-l)2tr(Q  VO  V) } 
a a e e 


and 


'2  2 

var  (o  ) = { tr  (Q  VQ  V) } . 

e (n-t>2 


III.  Results 

The  measure  of  sensitivity  of  the  variance  of  the  estimators  is 
defined  as  follows: 


variance  of  the  estimator  of  a being  considered 

— 1 

minimum  variance  in  class  of  quadratic,  unbiased  estimators 

2 

The  measure  used  does  not  involve  the  variance  of  the  estimator  of  a 

e 

being  considered  since  for  all  cases  investigated  the  variance  of  o 

e 

"2 

remained  stable  relative  to  the  variance  of  for  various  parameters, 

a priori  weights,  and  designs.  Note  that  y >1  and  that  the  decimal 

a “ 

portion  of  y times  100  is  the  percent  by  which  the  variance  of  the 

a 

2 

estimator  of  being  considered  exceeds  the  minimum  variance  for  a qua- 
dratic, unbiased  estimator.  Also,  y depends  on  the  variance  parameters 

2 a 

only  through  the  ratio  — . 

°e 

In  order  to  compare  results  obtained  across  different  sample  sizes 
and  different  number  of  levels  of  the  random  effect,  a measure  of  un- 
balancedness is  defined.  Using  an  extension  of  a measure  used  by  Low  [1976] , 
D 


% max  D - 


max 


and  D is  the 
max 


p £ *1  £ 

where  D - l (r  -r)2  - min  £ (r  -r)2 I,  r - y £ r = j , 

i-1  1 Li-1  J i-1  1 1 

maximum  value  D can  attain.  Designs  are  then  chosen  for  comparison  that 
have  approximately  the  same  % max  D.  Note  that  when  the  design  is  balanced 


6 


as  possible  then  % max  D = 0 and  that  when  the  design  is  as  unbalanced 

as  possible  then  % max  D = 1. 

The  study  took  many  design  configurations  into  consideration.  That 

°a  Wa 

is,  the  behavior  of  y for  values  of  — - , — , n,  i,  and  r.  (i=l,2, H) 

a i w l 

°e  e 

was  studied.  Designs  were  considered  where  more  degrees  of  freedom  were 

2 2 

available  to  estimate  o than  a ; that  is,  the  majority  of  the  r.  =1. 

a e l 

The  results  reported  here,  however,  are  only  for  the  case  when  there  are 

2 2 

more  degrees  of  freedom  to  estimate  o than  o . Attention  is  devoted  to 

e a 

this  type  of  design  for  two  reasons.  First,  this  seems  to  be  the  case 
in  the  majority  of  physical  situations.  Second,  there  is  a lack  of  con- 
sistency and  trends  when  more  degrees  of  freedom  are  available  to  estimate 
the  random  effect  variance  than  the  error  variance.  The  results  which 
are  tabulated  and  graphed  below  are  representative  of  the  many  design 

configurations  considered. 

[Insert  Tables  and  Figures] 

IV.  Conclusions 

When  considering  designs  for  which  there  are  more  degrees  of  freedom 
, 2 2 

to  estimate  o than  o and  for  which  the  unbalancedness  is  not  extreme, 
e a 

recommendations  can  be  made.  As  can  be  seen  from  tables  1 and  2 the 
MINQU  estimator  becomes  increasingly  sensitive  to  the  a priori  weights 
as  the  degree  of  unbalancedness  increases.  Also,  tables  1 and  2 as  well 
as  figures  1 and  2 illustrate  that  as  the  number  of  levels  of  the  random 
factor,  l,  increases  the  more  sensitive  the  MINWU  estimator  becomes  to 


the  ratio  of  the  weights. 


w 


w 


Again  looking  at  tables  1 and  2,  graphs 


1 and  2,  and  excepting  the  case  when  the  variance  ratio  is  smallest 

it  is  seen  that  the  MINQU  estimator  is  more  sensitive  to  the  a priori 

°a 

weights  which  underestimate  — - than  it  is  to  the  a priori  weights  which  in 

« l u 

ratio  are  greater  than  — . 

°e 


Incorporating  the  results  obtained  for  the  Henderson  analysis  of 
variance  type  estimator  (Table  3)  with  the  results  noted  above  for  the 
MINQUE  the  following  recommendations  can  be  made  for  the  one  way  random 


effect  model. 

a2 
a 

2 


variances. 


First,  if  the  experimenter  believes  that  the  ratio  of 
, is  less  than  one  then  the  Henderson  Method  I estimator 


should  be  used.  Second,  if  the  experimenter  believes  that  the  random 
effect  variance  is  greater  than  the  error  variance  then  the  MINQU  estima- 
tor should  be  used  incorporating  a priori  weights  which  in  ratio  you  feel 

°a 

does  not  underestimate  — . 

a 


e 
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Table  2.  MINQUE — Y for  various  parameters, 

d 

a priori  weights,  and  design  configurations 

2 


a 

2 


= 1.0 


w 

a . . . . 

— is  this  por- 

We  of 

tion  of  — 
a 

e 


•design 

*=3 

ri=4 

r2=6 

r3=8 

1=5 

ri=2  r4=5 
r2=3  r5=6 

r3=4 

«.=3 

ri=2 

r2=3 

r3=!3 

t=5 

V2  r4=2 
r2=2  r5=12 

r 3=2 

% max  D 
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3 
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w 
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— W 
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1.001 

tion 

r a 

2.00 

1.000 

1.000 

1.000 

1.000 

5.00 

1.000 

1.000 

1.001 

1.001 
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Captions  for  Figures 


Figure  1 

Y for  MINQUE  for  various  a priori  weights. 

u 

right  ordinate  indicate  the  constant  value  of  Y 

2L 

I estimator.  Design:  Jt  = 3,  r^  = 2,  r2  = 3,  r^ 


Figure  2 

y for  MINQUE  for  various  a priori  weights. 

a 

right  ordinate  indicate  the  constant  value  of  y 

a 

I estimator.  Design:  *.  = 5,  r^  = 2,  = 2,  r3 
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